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基于 计算 均匀 化 方法 的 梯度 多 孔 结构 优化 设计 


陈 炼 雄 , 肖 奥 宁 ,刘辉 


(武汉 大 学 土木 建筑 工程 学 院 工 程 力学 系 ,430072 武汉 ) 


Ho 要 :多 孔 结 构 凭借 轻 量 化 .高 比 强度 等 优越 性 能 ,在 现代 工业 领域 受到 越 来 越 多 的 关注 。 增 材 
制造 技术 的 快速 发 展 ,使 得 具有 复杂 构 型 的 多 孔 结 构 的 制备 成 为 可 能 ,进一步 促进 了 多 孔 结 构 设 计 

方法 的 进步 。 多 孔 结 5 构 设计 方法 主要 有 尺度 分 离 和 尺度 关联 两 种 途径 ， 本 研究 从 尺度 分 离 的 思想 

出 发 ,基于 计算 均匀 化 方法 ,针对 随机 (两 种 维 诺 型 ) 和 非 随机 ( 方 框 型 ) 梯 度 多 孔 结 构 , 发 展 了 Pus 
度 优 化 设计 方法 。 该 方法 主要 包括 3 个 步骤 :微观 尺度 微 结构 建 模 分 析 ( 涵 盖 微 结构 几何 生成 、 等 

效 物 理 量 的 计算 几何 参数 与 各 物理 量 之 间 映 射 关 系 的 建立 )、 宏 观 尺 度 结构 优化 和 全 尺度 多 孔 结 

构 几 何 重 构 。 该 方法 具 L 有 较 高 的 计算 效率 ,同时 保证 了 相 邻 微 结构 之 间 的 完美 连接 。 利 用 贴 体 三 

角形 网 格 划 分 ,建立 多 孔 结 构 全 尺度 有 限 元 模型 ,以 方 框 型 梯度 多 孔 结构 和 两 种 维 诺 型 梯度 多 孔 结 

wt | ,开展 了 多 孔 结 构 尺 寸 效 应 分 析 及 三 点 弯曲 梁 数 值 试验 ,验证 了 本 研究 方法 的 有 效 性 和 梯度 
多 孔 结 构 的 优越 性 ,研究 表明 多 孔 结 构 相 对 于 单 尺 度 结 构 具 有 和 较 好 的 承载 能 力 和 延展 性 

ean 梯度 多 孔 结构 ;多 尺度 优化 ;计算 均匀 化 ;尺度 分 离 ; 尺 寸 效应 
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Optimization design of graded cellular structure based on 
the homogenization approach 


CHEN Lianxiong, XIAO Aoyu, LIU Hui 


( Department of Engineering Mechanics ,School of Civil Engineering of Wuhan University ,430072 Wuhan, China) 


Abstract ; Cellular structures have attracted increasing attention in modern industry because of their excel- 
lent properties such as light weight and high specific strength. The rapid development of additive manufac- 
turing technology makes it possible to fabricate cellular structures with complex configurations , which fur- 
ther promotes the progress of cellular structure design methods. There are two main methods of cellular 
structure design :the scale-separated approach and the scale-related approach. Based on the idea of scale- 
separated approach and computational homogenization method , this paper develops a multi-scale optimiza- 
tion design method for stochastic (two kinds of Voronoi-shaped) and non-stochastic ( box-shaped) graded 
cellular structures. This method mainly containsthree steps,i. e. ,the modeling and analysis of microstruc- 
ture on the microscale ( including the generation of microstructural geometry , calculation of equivalent 


physical quantities , and establishment of mapping relationship between geometric parameters and physical 
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quantities) ,optimization of the macrostructure on the macroscale ,and reconstruction of the graded cellular 


structure on the full-scale. This proposed scale-separated method has high computational efficiency and en- 


sures the connectivity of the microstructure perfectly. The full-scale finite element model of the obtained 


cellular structure is established by using the triangle-fitting mesh division. The size effect analysis and the 


three-point bending beam numerical experiment of the cellular structures are carried out by taking the box- 


shaped graded cellular structure and two kinds of Voronoi-shaped graded cellular structure as examples, 


which verify the effectiveness of the proposed method and the superiority of the obtained graded cellular 


structure. The results show that the bearing capacity and the ductility of the graded cellular structure are 


better than the single-scale structure with little stiffness loss. 


Key words: graded cellular structure ; multiscale optimization ; homogenization ; scale-separated ;size effect 


随 着 科学 技术 的 不 断 发 展 ,人 类 对 新 材料 和 新 
结构 的 需求 越 来 越 高 ,高 性 能 、 低 能 耗 等 成 了 不 必 可 
少 的 条 件 。 在 工业 需求 的 驱动 下 ,结构 材料 一 体 化 
读 计 方法 研究 备 受 推崇 , 则 在 通过 多 尺度 设计 途径 
得 到 性 能 更 加 优越 的 多 尺度 结构 。 增 材 制造 技术 的 
快速 发 展 "1 ,进一步 解决 了 多 尺度 结构 生产 制造 难 
晤 5 雅 动 形成 多 尺度 结构 发 展 应 用 新 格局 。 
回 自 然 界 中 的 竹子 .骨头 等 多 孔 结 构 是 典型 的 多 
尽 庆 结构 ,具有 重量 轻 、 强 度 高 . 耐 冲击 等 优点 "1， 
得 到 了 工业 领域 工程 师 们 的 青睐 。 多 尺度 结构 通常 
惹 有 两 个 或 多 个 特征 尺度 :材料 尺度 称 为 微观 尺度 ， 
绪 多 尺度 称 为 宏观 尺度 。 两 个 尺度 相互 耦合 ,微观 
椒 料 的 分 布 会 影响 其 宏观 物理 性 能 ,宏观 结构 边界 
条 种 的 改变 又 会 引起 局 部 微观 材料 分 布 的 变化 。 实 
现 肖 微观 尺度 结构 和 材料 一 体 化 设计 是 极 具 挑 战 性 
的 课题 。 学 者 们 发 展 了 许多 具有 代表 性 的 多 尺度 方 
法 = 如 计算 均匀 化 方法 "1 多 尺度 有 限 元 方法 "中 和 
全 及 度 / 减 缩 尺度 方法 沾 。 这 些 方法 可 以 分 为 两 
类 , 即 尺度 分 离 途径 和 尺度 关联 途径 。 尺 度 分 离 途 
径 的 核心 思想 是 将 结构 划分 为 宏观 尺度 和 微观 尺 
度 , 基 于 计算 均匀 化 方法 或 代表 体 元 法 等 ,计算 微观 
尺度 的 宏观 等 效力 学 性 能 ,该 途径 具有 较 高 的 效率 ， 
但 是 需要 解决 微观 尺度 结构 连接 性 问题 ,并 开展 尺 
寸 效 应 分 析 。 尺 度 关联 途径 是 建立 宏观 和 微观 尺度 
物理 量 之 间 的 联系 ,并 进行 全 尺度 优化 设计 ,该 方法 
保证 了 微 结构 完美 连接 ,但 是 通常 计算 效率 较 低 , 难 
以 用 于 大 规模 多 尺度 结构 的 优化 设计 。 因 此 ,开展 
多 孔 结构 宏 微 观 一 体 化 设计 方法 研究 ,保证 计算 效 
率 和 微 结构 连接 性 ,具有 重要 的 意义 。 

按 微 结构 类 型 划分 ,多 孔 结构 可 分 为 非 随机 多 
孔 结构 "1 和 随机 多 孔 结构 ** ,其 中 非 随机 多 孔 
结构 以 方 框 型 X 型 等 基本 微 结构 构 型 及 其 组 合 构 


型 组 成 ,这 些 基 本 构 型 结构 形式 简单 ,几何 参数 较 
少 , 可 控 性 良好 ,但 同时 具有 较 强 的 各 向 异性 特征 ， 
导致 多 尺度 结构 对 边界 条 件 的 敏感 性 显著 增 大 。 而 
随机 多 孔 结 构 以 不 规则 多 边 型 微 结构 为 主 ,如 维 诺 
型 ,随机 微 结构 表现 出 较 好 的 各 向 同性 性 质 , 在 实际 
应 用 中 表现 出 优越 的 力学 行为 。 目 前 对 于 随机 微 结 
构 优化 设计 的 研究 并 不 多 , Wu 等 ” 基于 均匀 化 和 
参数 化 方法 ,提出 了 一 种 共 形 多 孔 结 构 优 化 设计 方 
法 ;Do 等 ”基于 均匀 化 途径 ,发 展 了 一 种 维 诺 随 机 
多 孔 结构 设计 方法 ;Lei 等 ”提出 了 一 种 参数 化 维 
诺 多 孔 结 构 设 计 方法 。 

本 工作 基于 计算 均匀 化 方法 ,针对 随机 和 非 随 
机 梯度 多 孔 结构 ,提出 了 相应 的 多 尺度 结构 优化 设 
计 方 法 。 该 方法 采用 了 尺度 分 离 思想 ,在 微观 尺度 
上 ,进行 微 结构 建 模 和 分 析 。 基 于 均匀 化 方法 ,计算 
微 结构 的 宏观 等 效 弹 性 张 量 ,建立 微 结构 物理 量 ( 宏 
观 等 效 弹 性 张 量 、 切 割 高 度 \、 平 均 长 细 比 ) 与 相对 密 
度 之 间 的 映射 关系 。 在 宏观 尺度 上 ,发 展 一 种 无 德 
罚 变 密度 拓扑 优化 方法 ,得 到 宏观 结构 相对 密度 分 
布 。 在 全 尺度 上 ,基于 宏观 相对 密度 分 布 和 微观 映 
射 关 系 ,计算 全 矿 度 结构 的 符号 距离 函数 或 切割 曲 
面 ,通过 切 制 符 号 距离 函数 ,得 到 全 尺度 多 和 孔 结构 几 
何 构 型 。 


1 微观 尺度 映射 关系 


本 节 将 介绍 非 随机 微 结构 和 随机 微 结构 物理 量 
和 相对 密度 之 间 映 射 关 系 的 建立 , 即 正则 化 近似 各 
向 同性 拉 梅 常数 、 正 则 化 弹性 矩阵 、 切 割 高 度 . 平均 
长 细 比 与 相对 密度 的 映射 关系 。 主 要 内 容 涉 及 微 结 
构 生 成 ,均匀 化 以 及 映射 关系 的 建立 等 。 其 中 非 随 
机 微 结构 以 方 框 型 微 结构 为 例 ,随机 微 结构 以 维 诺 
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型 微 结 构 为 例 。 
1.1 微 结构 生成 


如 图 1 所 示 , 随 机 维 诺 型 微 结构 生成 过 程 中 , 基 
于 网 格 生 成 算法 ,在 微 结构 域内 随机 生成 节点 ， 
形成 网 格 训 分 , 微 结构 域 边 界 节 点 需 满足 对 称 性 要 
求 ,不 断 和 迭代 并 收敛 形成 均匀 的 三 角形 节点 和 网 格 
剖 分 。 利 用 维 详 型 和 三 角形 网 格 放 分 在 空间 几何 的 
耦合 关系 ” ,得 到 维 诺 型 网 格 训 分 , 即 维 诺 型 微 结 
构 骨 架 。 基 于 快速 行进 方法 ” ,计算 得 到 微 结构 骨 
RWI FERRA, (x,y)。 通 过 水 平 集 函 数 $, = 
f. Gy) +to, 得 到 维 诺 型 微 结构 , 即 

,>0, 微 结构 实体 材料 


lo. =0, 微 结构 边界 a) 
> 中 ,<0, 微 结构 空 材料 


ju =1/2,t 是 微 结构 内 部 杆 件 的 宽度 。 


(e) 微 结构 


(b) 符号 距离 函数 


N 
d 图 1 随机 维 诺 型 微 结构 生成 过 程 
Ae Fig. 1 Generation process of the stochastic 
P Voronoi-shaped microstructure 
SEME 2 所 示 , 非 随机 方 框 型 微 结构 生成 过 程 中 ， 
KER GC EUNT. 
= e, »max(2|£| -1,2|m| - 1) (2) 


HE &( -1<e<1) 和 1( -1<m<1) 是 微 结构 域 局 
MUI. 

用 不 同 高 度 平面 切割 水 平 集 函 数 ,将 得 到 不 同 
方 框 型 微 结构 。 由 于 水 平 集 函 数 的 取 值 范围 为 
[ -1,1] , 当 切 割 高 度 小 于 -1 时 ,将 得 到 实体 方 框 
型 微 结构 ; 当 切 割 高 度 大 于 1 时 ,将 得 到 空白 方 框 型 
微 结构 ;图 2(e) 中 所 示 为 切割 高 度 为 0.6 时 对 应 的 


方 框 型 微 结构 。 
n. 2 
E E -0.5 
: 14 EN 1 
70.0 
ST ki 
(a) 微 结构 域 (b) 水 平 集 函 数 (c) 微 结构 


图 2 非 随 机 方 型 微 结构 生成 过 程 


Fig.2 Generation process of the non-stochastic 


"o 
陈 炼 峻 ， 等 ， 基 于 计算 均匀 化 方法 的 梯度 多 孔 结 构 优化 设 训 inaXiv 合 作 期 刊 1057 


box-shaped microstructure 
1.2. 微 结构 均匀 化 


根据 计算 均匀 化 理论 ,周期 性 微 结构 的 宏观 


H Olij) (ij) OCKL) (kL) 
E, = | Y | Í, E. ( £p Hom ES ) ( £, 一 Cr ) dY 


(3) 
Hop: LY, ifte Mg SE, (x,y) 是 局 部 变 刚 度 
弹性 张 量 ;eM 是 局 部 变化 的 应 变 场 ;es” 是 已 知 的 
宏观 应 变 场 。 对 于 二 维 问题 , e ARN ES = 
(1,0,01, 227? = {0,1;0 关 和 eg) = (0,0,1]*, 
根据 Voigt 记号 , 微 结 构 等 效 弹性 张 量 Eji, 可 以 
写成 矩阵 的 形式 
Di D D 
D' =| Dì) Do XLI =D} (4) 
Ds Da D 
表 1 给 出 了 一 个 维 谤 型 微 结 构 和 方 框 型 微 结 构 
等 效 弹性 矩阵 ,其 中 微 结构 实体 材料 和 空 材料 的 无 
A 34 du dont X Ap 5] DJ (Alu 2.25 4) 和 
(A, 710 A ad ;Hvoid = 10 ue 从 表 1 可 以 看 
出 , 方 框 型 微 结 构 等 效 弹性 矩阵 只 有 3 个 参数 , 维 诺 
型 微 结构 等 效 弹性 矩阵 有 6 个 参数 ,其 中 有 2 个 参 
数值 很 小 。 进 一 步 地 , 当 维 诺 型 微 结 构 内 部 杆 件数 
量 足 人 够 多 时 , 微 结构 将 趋 于 各 向 同性 ,因此 有 必要 对 
维 诺 型 微 结构 各 癌 异 性 弹性 进行 近似 各 向 同性 化 
处 理 。 
表 1 均匀 化 结果 


Tab.1 Results of the homogenization 


p! 
1.4380 0.8601 0.007 2 
[ose 1.4397 ni 
0.0072  -0.0021 0.3056 


(CC 
1.9993 0.1122 0 
|o: 1.9993 0 | 
0 0 0.0544 


近似 各 向 同性 理论 ”是 将 各 向 异性 弹性 张 量 
空间 通过 正 交 算 子 , 映 射 到 各 向 同性 张 量 子 空间 。 
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基于 此 ,近似 各 向 同性 拉 梅 常数 XA“ 和 jy 计算 如 下 。 


iso 1 
A” =y (Pi + DÀ, -ADE +6D5) 


(5) 
iso 1 
u -g (Di +D}, «ADS, -2Di) 
AERIS BE EUR i [I] ERR ED" 
E 124^ je 0 
人 Di - Ae Qj se 124^ 0 (6) 
0 0 u™ 


为 了 量化 近似 各 向 同性 化 处 理 的 误差 ,各 向 异 
性 指数 六 定义 为 


2 


ID 


L-5h- (7) 
T |." |? 

Hep lp l3 D" 下 表达 式 为 

b mm 


2 二 4 (33s +12 Quy a 8 jj 


| D^" 

3 , , (8) 

ID" = pa (DD wd (DH)? 

加 结合 式 (5) (7) 和 (8) ,得 到 表 1 中 维 诺 型 微 结 

梅 芙 效 弹性 矩阵 各 向 同性 指数 1, =0.011, 接 近 于 0， 

说 明 近似 各 向 同性 化 处 理 是 较为 合理 的 。 

T 

13 微 结 构 映射 关系 建立 

忌 计算 微 结构 相对 密度 时 ,采用 了 固定 背景 网 格 。 

如 图 3 所 示 , 微 结构 中 有 3 种 网 格 单元 :实体 单元 、 

空 丰 单元 和 混合 单元 。 微 观 结构 的 相对 密度 p,, 的 计 
下 。 


c f Hcg)ay pA 
Vm i 
Pn = y D» (9) 
Q dV m | i=l 


Vm 
其 中 :H(x) 是 Heaviside PR ; V,, XAR WON Aa 3] R; 
NN. 表 示 单 元 数量 ;p 是 微观 结构 中 第 i 个 单元 的 材 
料 体积 。 


图 3 背景 网 格 示意 图 
Fig.3 Illustration of the background grid 


pw 的 表达 式 为 


ChinaXiv 合 作 期 刊 第 39 郑 


JiC£ D (w, 


p, = | HCA = È Hong) 

l (10) 

Hp J.B i MAAIE n EGXB PEE, dS W 

积分 点 的 坐标 向 量 ;w, 表示 积分 权重 因子 ;VV, 为 第 i 

个 单元 域 ;NN, 为 高 斯 积分 点 个 数 。 

在 维 诺 型 微 结构 物理 量 与 相对 密度 映射 关系 建 

立 过 程 中 ,建立 了 3 种 映射 关系 。 通 过 改变 切割 高 

度 ,计算 对 应 微 结构 的 相对 密度 ,得 到 切割 高 度 h 与 

相对 密度 之 p, 间 离散 的 映射 关系 ,并 采用 样 条 插值 
得 到 连续 的 映射 关系 ,如 图 4(a) 所 示 。 

0.10 


0.08 


0.06 


度 A 


Ir 
T E 
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20 

a 拟 合 数据 
拟 合 曲线 
15 E 验证 数据 


1, Jt 


g= 


(6) 平均 长 细 比 


1.0 
4 AG UE 
一 一 4 拟 合 曲线 
0.8 六 4 验证 数据 R 
z a LM $ 
E --- 4/ 拟 合 曲线 r 
ge 0.6 全 /验证 数据 ,所 
之 
E 0.4 
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(c) 正则 化 近似 各 向 同性 拉 梅 常数 
图 4 维 诺 型 微 结 构 物 理 量 与 相对 密度 映射 关系 


Fig.4 Mapping relationship between relative density and 


physical quantity of Voronoi-shaped microstructure 
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同样 地 ,采用 指数 函数 拟 合 内 部 杆 件 平均 长 细 

比 a=7/it 与 相对 密度 p,, 的 关系 ,如 图 4(b) 所 示 ， 
这 里 1,,, 为 杆 件 平均 长 度 , 且 

o = ae" +ce 

其 中 :a,b,c 和 是 拟 合 参数 ,如 表 2 所 示 。 
R2 MES a,b,c Hd 


Tab.2 Fitting parameters a,b,c and d 


dpm 


(11) 


a b c d 


43.8 - 13. 46 10. 15 -2.107 


此 外 ,采用 多 项 式 拟 合 正 则 化 近似 各 向 同性 拉 


梅 d 数 ( À norm m À AA solid sMnorm = uu us) 与 相 对 密 
度 p, 之 间 的 关系 ,如 图 4(e) 所 示 。 


Aum 7 Pip» PPa PPh 十 下 po 十 站 
| 4 3 2 (12) 
Unom = Pis + Pops + spa + Dips + Ds 

HEP p^ (i 21-5) DA SC, WME 3 所 示 。 

O e Jr ftt RU BL AG CRESCE vr ERE rh, 如 图 
5( 所 示 , 采 用 多 项 式 拟 合 方 型 微 结构 正则 化 等 效 
TEPEXEPED; (on) - Dj (p,)/D, (1) (六 =11,12,33)， 
EH D? (1) 是 实体 材料 (p。 =1) 单 元 对 应 的 弹性 
答 阵 。 


1.0 

dig o 万, 拟 合 数据 
i» D; 拟 合 曲线 
0.8 " Dj ru jE 
i 4 Di e RR 

-3 0.6| 77 - DHL HR 

c 人 A 站, 验证 数据 

A + 及, 拟 合 数据 

aj 04 Di MA d£ 

E x D 验证 数据 


> 
to 


0.0 


(a) 正则 化 等 效 阐 性 矩阵 


图 5 
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AH 11 4 11 3 11 2 11 
Du(p,) =a p, +D Pm 十 C p, *d Pn 
DH 12 4 12 3 位 5 12 
Dolpa) -a p» +D Pm tC p, td Pm (13) 
MH 33 4 33 3 33 2 33 
D; (Pan) = a Pm + b Pm + C Pn 十 d Pm 
其 中 ju" 2 Nd ,用 是 拟 合 参数 , 如 表 4 所 示 。 
X3 Hiep T (i =1 ~5) 
Tab.3 Fitting parameters p? and pr (i=1~5) 
i=l i-2 i-3 i-4 i-5 
p 2.912 -4.545 1.189 1.047 0. 009 
DA 2.192 -2.276 2.158 -0.294 0.015 
表 4 HUS Vb! cF d’ 
Tab.4 Fitting parameters a ,b? ,c/ and d? 
al bï c? d 
y=11 1.5270 -1.9830 1.1020 0.350 4 
ij -12 0.4695 0.4783 — 0.009 9 0.066 4 
ij 233 1.9340 —1.2610 0.3419 —0.0134 


如 图 5(b) 所 示 , 方 框 型 微 结构 切割 高 度 有 与 相 
对 密度 p,, 之 间 的 连续 映射 关系 通过 样 条 插值 得 到 。 


0.5 
0.0 
20:5 


0.8 


(b) 切割 高 度 


量 与 相对 密度 映射 关系 


Fig.5 Mapping relationship between relative density and physical quantity of box-shaped microstructure 


2 宏观 优化 列 式 


本 研究 提出 了 一 种 无 惩罚 密度 优化 方法 ,考虑 
最 小 柔顺 性 优化 问题 ,优化 列 式 为 
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min:J(p) = F"U 
KU =F 
V-V<0 


P min <p SP max 


OPT 


(14) 
s.t. ， 
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其 中 :7 是 目标 函数 ;p 是 相对 密度 ;是 外 力 荷载 ;U 
是 位 移 向 量 ;是 刚度 矩阵 ; V = [pdQ 是 结构 体积 


分 数 ;V 是 结构 体积 分 数 上 限 ;p,, 和 pi 分 别 为 相对 
密度 上 下 限 。 
在 优化 过 程 中 ,基于 有 限 元 方法 对 结构 进行 离 
散 计算 ,采用 OC 准则 求解 优化 问题 。 
目标 函数 对 相对 密度 的 灵敏 度 为 
aJ _ _yrðK 


15 
op, T Vd 
eN K, 
对 于 维 诺 型 多 孔 结构 ,3o 计算 如 下 。 

Ne 2 

oK € T n oD™ 
—-Yyc([m--naQ)c, (16 
0p, 2 (J, ð p. o) ue, 


Hah: C, 是 扩充 矩阵 ,用 以 将 单元 刚度 集成 到 总 


Hl B 是 应 变 箱 阵 ;人 是 单元 域 ;一 计 算 如 下 。 


OA +2 ðu QA 0 
- 9 p, 9 p, ð p, 
D” - 3 AP ð A " 2 9 iso T 
9 p, à p, ð pe à p. 
G | 0 0 ðu 
N (17) 
Hap 
| 0A E 9A norm 
óp, = solid óp, , 
ou^ ES N OL norm 1 8 
= ap Msolid ap ( ) 
dA 
m =A aia (4 pip; +3 pip, *2 pip, * pi) 
: (19) 
ÕM norm E 4 u 3 3 u2 2 m m 
Fr = aia (4 Pipe +3 pop, *2 psp, * pi) 


对 于 方 杠 型 多 孔 结构 ,25 计算 如 下 。 


Y c(| B ¿P Bao)c, (20) 


ð p, ES! 
其 中 

aD, Dij, " 

op, Pp 

H [gD* gop! 

oD E 12 11 0 (21) 

op op, Op 
aD} 
0 33 
op, 


| 
学 学 报 adul TI 第 39 着 
óDi ôD u 
à 三 Dj ( 1 ) 
Pe op。 
api. aD" 
e LIED) (22) 
op, — Op, 
am anu 
à = D33 ( 1 ) 
p. 9p. 
aD =4a" pi +3 "p? +2 cp, +d!!! 
op, 
aD = 12 3 12 2 12 12 23 
——- =4a “p, *3b"p, *2c p, +d (23) 
op, 
aD” 
P 2 =4a p +3b” p? +2c”p, id 
p. 


对 于 方 框 型 和 维 诺 型 多 孔 结构 ,目标 函数 对 体 
积 约束 的 灵敏 度 为 


oV _ H 
n. = han = y, (24) 
其 中 : 亿 表 示 单 元 体积 4 数 。 


3 ”全 尺度 多 孔 结 构 几 何 重 构 


本 小 节 以 图 6(a) 所 示 的 悬臂 梁 为 例 , 尺 寸 为 
W/L=17/2 ,介绍 了 维 诺 型 梯度 多 孔 结 构 及 方 框 型 梯 
度 多 孔 结 构 全 尺度 几何 构建 过 程 ,其 中 维 诺 型 梯度 
EILE I DORAT 2 种 方式 进行 构建 。 方 框 
型 多 孔 结构 和 工 型 维 诺 多 孔 结 构 是 利用 切割 高 度 与 
相对 密度 映射 关系 重 构 的 , 工 型 维 诺 多 孔 结构 是 基 
于 杆 件 平均 长 细 比 与 相对 密度 的 映射 关系 重 构 的 。 

如 图 6(b) 和 图 6(c) 所 示 , 对 于 方 框 型 梯度 多 
孔 结构 和 1 了 型 维 诺 梯 度 多 孔 结 构 , 得 到 宏观 尺度 优 
化 相对 密度 分 布 后 ,基于 微观 尺度 切割 高 度 和 相对 
密度 的 映射 关系 ,得 到 切割 高 度 分 布 ,构建 切割 高 度 
曲面 切割 结构 域 水 平 集 函 数 ,得 到 梯度 多 孔 结 构 
构 型 。 

如 图 6(d) 所 示 , 对 于 开 型 维 诺 型 梯度 多 孔 结 
构 ,得 到 宏观 尺度 优化 相对 密度 分 布 后 ,基于 微观 尺 
度 内 部 杆 件 平均 长 细 比 与 相对 密度 的 映射 关系 ,得 
到 平均 长 细 比 分 布 。 基 于 平均 长 细 比 分 布 ,借助 网 
格 生 成 算法 ,得 到 J 型 维 庄 多 孔 结构 骨架 。 基 于 快速 行 
进 方法 ,计算 骨架 的 符号 距离 函数 ,切割 基于 骨架 符号 
距离 函数 构建 的 水 平 集 函 数 , 得 到 梯度 多 孔 结 构 。 
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(c) I 型 维 诺 梯度 多 孔 结 构 


00130v1 


4 数值 算 例 
= 


CX p gage tg c e o H, T RR REA 
析 方 法 ,开展 了 尺寸 效应 研究 和 三 点 弯曲 梁 数 值 斌 
验 ` 验 证 本 研究 所 发 展 方法 的 正确 性 ,探究 基于 该 方 
法 移 化 设计 的 方 框 型 和 两 种 维 诺 型 梯度 多 孔 结构 性 
能 后 优越 性 。 其 中 尺寸 效应 分 析 基 于 如 图 6(a) 所 
示 辆 悬臂 梁 模 型 ,三 点 弯曲 梁 数值 试验 基于 如 图 7 
所 东 的 简 支 梁 模型 ,尺寸 W/L=1/3。 考 虑 平面 应 变 


门 是 所 有 物理 量 都 是 无 量 纲 的 。 


F 


L 


图 7 三 点 弯曲 梁 数 值 试验 模型 


Fig.7 Model of the three-point bending numerical experiment 


L, 


L, 


4.1 尺寸 效应 分 析 


本 小 节 通 过 对 比方 框 型 多 孔 结构 ( 微 结构 数量 
分 别 为 60 x30,30 x 15,20 x 10) , 工 型 维 诺 多 孔 结 
构 ( 胞 元 数量 分 别 为 2 000 1 000.500) ,] 型 维 诺 多 孔 
结构 (内 部 杆 件 宽度 分 别 为 0.032 .0.043 6,0. 055 4) ， 
探究 所 发 展 的 多 尺度 方法 在 优化 多 孔 结 构 时 的 尺寸 


ea 
陈 炼 峻 ， 等 ， 基 于 计算 均匀 化 方法 的 梯度 多 孔 结 构 优化 设 训 inaXiv 合 作 期 刊 1061 


(b) 方 框 型 梯度 多 孔 结构 


(d) 工 型 维 诺 梯度 多 孔 结构 


图 6 梯度 多 孔 结构 重 构 


Fig.6 Reconstruction of the graded cellular structures 


效应 。3 种 多 孔 结 构 的 密度 上 下 限 均 为 p,,, =0.8 和 
p; =0. 20 

如 图 8 所 示 , 借助 Abaqus 软件 ,基于 全 尺度 有 
限 元 分 析 , 与 宏观 优化 结果 对 比 ,验证 方法 的 正确 
性 。 在 宏观 尺度 优化 过 程 中 , 方 框 型 多 和 孔 结构 、I 型 
和 本 型 维 诺 多 孔 结 构 的 柔 度 值 分 别 为 13.31、13.23、 
13. 19。 全 尺度 有 限 元 计算 得 到 的 3 种 方 框 型 多 孔 
结构 ( 微 结构 数量 分 别 为 60 x 30,30 x 15,20 x 10) 
承受 单位 外 荷载 时 的 最 大 位 移 分 别 为 12.33、13. 08、 
13.67,3 种 工 型 维 诺 多 孔 结 构 ( 微 结构 数量 分 别 为 
2 000 1 000.500) 承 受 单 位 外 荷载 时 的 最 大 位 移 分 
别 为 12. 53,12. 91,13. 48,3 $h I Æ AE £ FLH r 
(内 部 杆 件 宽 度 分 别 为 0. 032,0. 043 6,0. 055 4) 承 
受 单位 外 荷载 时 的 最 大 位 移 分 别 为 12. 19,12. 73, 
13.24。 表 5 给 出 了 有 限 元 结果 与 宏观 优化 结果 误 
差 ,其 中 “-” 表 示 有 限 元 结果 比 优 化 结果 小 ,“+” 
表示 有 限 元 结果 比 优化 结果 大 。 可 以 看 出 ,全 尺度 
有 限 元 计算 的 结果 与 宏观 尺度 优化 结果 比较 接近 ， 
表明 本 研究 所 提出 的 多 孔 结 构 优 化 的 设计 方法 是 正 
确 性 与 有 效 性 。 

通过 对 比方 框 型 多 孔 结 构 不 同 微 结 构 数 量 、[ 
型 维 诺 多 孔 结 构 不 同 微 结 构 数 量 及 开 型 维 诺 多 孔 结 
构 不 同 内 部 杆 件 宽度 全 尺度 有 限 元 计算 过 程 中 结构 
承受 单位 外 荷载 时 的 位 移 云 图 ,可 以 看 出 ,位 移 分 布 
近乎 相同 。 通 过 对 比 最 大 位 移 值 ,可 以 看 出 , 随 着 多 
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孔 结 构 微 结构 数量 增多 或 内 部 杆 件 宽度 变 小 ,多 和 孔 

结构 更 加 密实 ,结构 刚度 更 大 ,对 应 的 结构 最 大 位 移 

也 逐渐 减 小 ,但 位 移 值 相 差 不 大 ,表明 基于 本 研究 方 

法 得 到 的 3 种 梯度 多 孔 结 构 均 对 尺寸 敏感 性 较 低 。 
表 5 全 尺度 有 限 元 尺寸 效应 分 析 误 差 


Tab.5 Error of full-scale finite element size effect analysis 


误差 方 框 型 /% HEWA I)/96 维 诺 型 ( )/96 
60 x 30/2 000/0. 032 -7.36 -5.29 -7.58 
30 x 15/1 000/0.0436 -1.76 -2.41 -3.49 
20 x 10/500/0. 055 4 +2.70 +1.89 +0.38 


ze 
— HE : 


i 


T 
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co 

T 
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e 

CN 

T 

CN 

Es (c) FHER 20x 10 
N 

2. MENT 
^c 三 点 弯曲 梁 数 值 试验 


CO inp 9 所 示 , 将 3 种 多 孔 结 构 与 基于 参 化 水 平 
集 咯 法 优化 得 到 的 单 尺 度 结构 作对 比 ,通过 三 点 弯 
数值 试验 ,探究 多 孔 结构 性 能 的 优越 性 。 在 全 

尺 厦 有 限 元 模拟 过 程 中 ,采用 了 超 弹 性 材料 本 构 关 


A 


尺度 结构 


工 型 维 诺 多 尺度 结构 


方 框 型 多 尺度 结构 
(a) 加 载 前 后 的 结构 构 型 


图 9 


Fig.9 Results of three-point bending numerical experiment 
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(£f) I Ei 500 
图 8 不 同 梯度 多 孔 结构 全 尺度 有 限 元 尺寸 效应 分 析 


Fig.8  Full-scale finite element size effect analysis of different cellular structures 
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需要 指出 的 是 , 当 相 对 密度 下 限 取 更 小 的 值 时 ， 
宏观 密度 分 布 中 相对 密度 小 的 区 域 ,如 结构 右上 和 角 
和 右 下 角 部 分 , 重 构 得 到 的 多 孔 结构 较为 稀疏 。 当 
相对 密度 下 限 趋 于 0 时 ,宏观 密度 分 布 中 相对 密度 
接近 0 的 区 域 , 重 构 得 到 的 多 孔 结构 几乎 为 空 材料 ， 
导致 局 部 结构 不 稳定 ,影响 整体 结构 性 能 。 因 此 本 
研究 相对 密度 下 限 取 为 0.2。 

不 可 否认 的 是 ,由 于 本 研究 密度 下 限 取 0.2, 结 构 
右上 角 和 右 下 角 部 分 对 刚度 贡献 不 大 的 区 域 , 仍 会 出 
现在 结果 中 ,需要 进一步 改进 ,以 期 提高 结构 性 能 。 


(e) IL AE: 0.032 


(h) IUE AES 0.043 6 


AHE 0.055 4 


系 , 考 虑 了 几何 非 线性 。4 种 结构 体积 分 数 上 限 均 
为 0.4, 方 框 型 多 孔 结 构 微 结构 数量 为 96 x30, 工 型 
维 诺 多 孔 结 构 的 微 结 构 数 量 为 2 000, 开 型 维 诺 多 孔 
结构 内 部 杆 件 宽度 为 0.24。3 种 多 孔 结 构 的 密度 上 
下 限 均 为 p， =1.0 flip, =0.2。 


一 一 单 尺度 结构 
M EE ERA 
— 型 维 诺 多 筷 结构 
Yaris 方 框 型 多 孔 结 构 e 
et 
5 Ce zx 
i 
"di 


位 移 /mm 


(b) 试验 过 程 力 位 移 曲线 


三 点 弯曲 梁 数 值 试验 结 
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4 种 结构 在 加 载 前 后 的 几何 构 型 如 图 9(a) 所 
示 , 可 以 看 出 , 单 尺度 结构 容易 发 生 整 体 失 稳 ,多 孔 
结构 内 部 发 生 局 部 失 稳 ,而 不 容易 发 生 整 体 失 稳 。 
图 9(b) 为 试验 过 程 的 力 位 移 曲线 ,可 以 看 出 ,多 和 孔 
结构 的 承载 力 和 延展 性 较 好 于 单 太 度 结构 , 单 矿 度 
结构 刚度 略 优 于 多 孔 结 构 , 表 明 多 孔 结构 在 刚度 损 
失 不 大 的 前 提 下 ,可 以 提高 结构 承载 力 和 延展 性 ,在 
承载 .抵抗 冲击 \ 吸 能 等 方面 具有 较 大 的 应 用 前 景 。 


5 结 论 


本 研究 基于 计算 均匀 化 方法 ,发 展 了 一 种 尺度 
分 离 的 梯度 多 孔 结构 多 尺度 优化 设计 方法 。 在 微观 
尺度 ,生成 了 微 结构 几何 构 型 ,基于 均匀 化 理论 计算 
了 微观 结构 的 宏观 等 效 弹性 性 质 ,建立 了 微 结构 物 
理 荀 (宏观 等 效 弹性 张 量 、 切 割 高 度 及 平均 长 细 比 ) 
与 相 对 密度 之 间 的 映射 关系 ;在 宏观 尺度 ,基于 无 征 
罚 的 变 密度 拓扑 优化 方法 ,得 到 了 宏观 结构 相对 窗 
诬 仅 布 ;在 全 尺度 几何 构建 过 程 中 ,通过 切割 结构 域 
捕 湾 平 集 函数 ,得 到 梯度 多 孔 结构 。 该 方法 具有 较 
高 的 计算 效率 ,同时 保证 了 相 邻 微 结构 较 好 连接 。 
数值 算 例 以 非 随机 方 框 型 多 孔 结构 和 随机 维 诺 型 多 
也 结构 ( 工 、 工 ) 为 例 , 与 传统 参数 化 水 平 集 方法 优 
化 得 到 的 单 尺度 结构 进行 对 比 , 通 过 全 尺度 非 线性 
有 和 限 元 分 析 , 探 究 了 多 孔 结构 的 尺寸 效应 和 三 点 弯 
晶 染 数值 试验 ,验证 了 本 研究 方法 的 有 效 性 及 梯度 
多 入 结构 性 能 的 优越 性 。 结 果 表明 ,梯度 多 孔 结构 
对 及 十 效应 敏感 性 较 低 ,承载 力 和 延展 性 相对 传统 
单 翁 度 结构 有 了 较 大 的 提升 。 
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